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For a weakly interacting Bose-Einstein condensate in a double well, an appropriate time-dependent 
modulation of the trapping potential counter-acts the "self-trapping" effects of the interactions, 
thereby allowing tunneling between the wells. It is demonstrated numerically that this modulation 
can be employed for transferring the condensate from one well to the other in a controlled way. 
Moreover it allows the production of mesoscopic entangled states on short time scales. 
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I. INTRODUCTION 

Bose-Einstein condensates in double well potentials 
display a wide variety of important phenomena, rang- 
ing from interference of two condensates released from a 
double well trap and the relative phase of two conden- 
sates over the Josephson effect ^^^^ to paramet- 
ric resonance 0- In this paper, we demonstrate that 
the tunneling effect between the two condensates can 
systematically be manipulated by an appropriate time- 
dependent variation of the double well potential. 

For a single particle in a double well, it has been shown 
that the tunneling effect can be coherently destroyed by 
an external oscillating force A similar effect has also 
been predicted for condensates 0- Here, we pursue the 
opposite direction: We demonstrate numerically that un- 
der the influence of a suitably designed force tunneling of 
condensates can take place even under conditions where 
it normally is inhibited by self-trapping. 

Self-trapping has only recently been observed experi- 
mentally for condensates in a double well potential [lol |: 
the ability to control tunneling in such systems will en- 
able the experimental realization of new phenomena. 
As a particular application of our technique, we pro- 
duce numerically highly entangled mesoscopic superposi- 
tions (an entangled state cannot be written as a product 
state). While possible applications of mesoscopic entan- 
gled states include interferometry, frequency standards 
and quantum information processing, the realization of 
mesoscopic superpositions of condensates is still a chal- 
lenge of fundamental research. So far, superpositions 
of single atom states [ll] |. four wave mixing of matter 
waves 

m and even m~e complicated quantum states 
showing multi-particle entanglement |l3j have been re- 
alized experimentally. There have been several sug- 
gestions to produce mesoscopic entanglement with con- 
densates by manipulating the interaction between the 
particles, o r by controlling the dynamics of the sys- 
tem [H [ll El [Hill El . The method presented here 
is quite simple, and robust against small uncertainties in 
both the preparation of the initial state and specific de- 
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tails of the force. It also might lend itself to a systematic 
adaption of highly successful methods developed in the 
field of coherent control of molecular dynamics poL I2H 
to condensates. 

Our material is organized as follows: First we intro- 
duce the A^-particle model used to describe the dynamics, 
and transform the A^-particlc Schrodinger equation into 
an equivalent set of equations. These equations reveal 
the guiding principle how to design the required time- 
dependent forces (a discussion of the validity of the model 
can be found in the appendix) . We then demonstrate nu- 
merically how an efficient population transfer from one 
well to the other can be achieved, and outline the strat- 
egy of creating mesoscopic entanglement. 

II. MODEL 

For a matter-of-principle discussion, we consider the 
iV-particle Hamiltonian 

+ hfif{t)(blk~blk) , (1) 

where the Bose operators Up annihilate (create) one par- 
ticle in the j-th well (j — 1,2). In addition, fi denotes 
the single-particle tunneling frequency, fin is the on-site 
interaction energy, and hfj,f{t) specifies an externally ap- 
plied potential difference between the two wells. Without 
this time-dependent modulation, this Hamiltonian is well 
established for describing the tunneling dynamics of a 
condensate in a double well potential within a two-mode 
approximation |2l . The time-dependent force has to 
be chosen such that this approximation is not violated. 

Considering two harmonic wells with angular fre- 
quency uj and oscillator length ^osc = {2mLu) sep- 
arated by a distance of 2c?, conservative estimates |23l 

H yield n ~ wy27^d/£oscexp(-0.5dV^osc) and k ~ 
ha / {Ay/TTmi^^^) , where m is the atomic mass and a 
the s-wave scattering length. Taking ^'^Na atoms with 
a — 2.75 nm, we obtain Nk ~ 1.651 for = 1000 par- 
ticles by choosing ^osc — 31.78/zm and d = 3.2£osc, say, 
with u ~ 1.362 Hz, 17 ~ 0.02078 Hz. Experimentally, 
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time-scales shorter by two orders of magnitude can eas- 
ily be obtained by either choosing just N = 100 or by 
reducing the interaction with a Feshbach resonance by a 
factor of 10. In both cases one would have: uj ~ 136.2 Hz 
and n ~ 2.078 Hz with 4sc ^ 3.178//m and d = 3.24sc- 
Without interaction between the particles {k — 0), 
and without any applied force (nfit) = 0), an initial 
state with all particles trapped in one well simply leads 
to Rabi-like oscillations of the entire population between 
both wells: {b\b^)/N = (1 ± cos(m))/2. For describing 
the evolution of the interacting system, we use the Fock 
basis = \N — i/, v) with v = ... TV, so that the la- 
bel V refers to a state with N — v particles in well 1, and 
V particles in well 2. The Hamiltonian (^) then becomes 
a sum of two {N + 1) x {N + l)-matrices, 

H = Hoit)+Hi. (2) 

The non-diagonal matrix Hi describes the tunneling. 



{i^\Hi\n) ^ ~ iS,,n-i + S.-i,n) —VN^V^TTl , (3) 

whereas the diagonal matrix Hq includes both the inter- 
action between the particles and the applied force, 

{iy\Ho{t)\n) - ,5,,„[(7V-z/)(7V-j/-l) + i/(i/-l)]fi,K 
+ S.^niN -2u)h^lf{t) . (4) 

To solve the Schrodinger equation 



we employ the ansatz 
{iy\ip{t)) = a^{t) exp 



(5) 



(6) 



A few lines of simple algebra then lead to 

iha^{t) = (iy\Hi\i^+l)K{t)at,+i{t) 

+ (i.|iJi|j.-l)/i,_i(i)*a,_i(t) 



(7) 



where we have set a_i(t) = ajv+i(^) = , and introduced 
the phase factors 
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FIG. 1: Maximum fraction of particles found in well 2, after 
initially all particles had been prepared in well 1, if no force 
is applied. Solid line: Prediction from Eq. for A'^ ^ 1. 
Dashed line: data obtained from the solution of the non-linear 
Schrodinger equation. There is a self-trapping transition for 

NK = n IMlil. 



condensates 0. Self-trapping can also be understood 
within the A'^-particle approach (c/. Ref. jl^l): Here, we 
start with a perfectly trapped situation, in which a^(0) = 
6^^Q. In the spirit of time-dependent perturbation theory, 
we then have to zeroth order ai,{t) = Si,^o for all times. 
In first order perturbation theory, we obtain ao{t) ~ 1, 
02 (t) ~ 0, a3(t) ~ 0, . . . , and 
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ai{t)~i-^/N exp[-i2(iV- l)Kt']di' 
2 Jo 



(9) 



Therefore, for not too large ^/[{N — 1)k], the fraction of 
particles in well 1 evolves in time as 



N 



1 - 



2{N -1)k 



\[{N - l)Ktf 



(10) 



so that the transfer of population is suppressed. 
Surprisingly, this perturbative reasoning qualita- 
tively describes self-trapping already for comparatively 
large ^/[{N — 1)k]; for Nk > 1.5 O even the quantitative 
agreement is good (see Fig.^. 



hv (t) = exp <^ i 



2{N - 1 - 2v)Kt + 2/i / f(t') dt' 



(8) 

This set of differential equations is mathematically equiv- 
alent to the iV-particle Schrodinger equation governed by 
the Hamiltonian It will enable us to design force 
functions /i/(i) to achieve the desired control of tunnel- 
ing. These predictions will then be tested numerically. 

If initially all particles are stored in one well, and 
no force is applied (/i/(t) = 0), tunneling is inhibited 
for Nk/U, > 1. This self-trapping effect, which origi- 



nally has been established on the mean- field level 0, |2E 
, has recently been observed experimentally for 
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III. TRANSFER OF BOSE-EINSTEIN 
CONDENSATES 

We now consider the question how to transfer the con- 
densate from one well to the other even under conditions 
of self-trapping. If there is no interaction, and all par- 
ticles are confined in well 1 at t = to, the ai,{t) which 
are essentially different from zero are centered around 
lydt) — N {1 — cos[n,{t — io)]) /2, as a consequence of the 
unperturbed Rabi oscillation. To reproduce this behavior 
in the interacting case, even in the self-trapping regime, 
the phase factor hi,^ has to remain constant at least for 
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small time differences At. From Eq. this requirement 
leads to the condition 

l~exp{i{2NKCOs[n{t-tQ)] + 2^f{t))At} . (11) 

To transfer a condensate initially prepared in well 1, we 
therefore design the force 

{-^J■Q ■ t <to 
cos[il{t - to)] : to <t <to + tt/O , 
fio ■ t >to + 

(12) 

with fii sa ~Nk. Fortunately, the resulting transfer is 
reasonably insensitive with respect to the precise value 
of III, as the system still performs Rabi-like oscillations 
for small but finite interactions. As in the experiments 
reported in Ref. 0|, the role of the constant offset 
is to prevent tunneling before and after the controlled 
evolution of the condensate. This offset is an auxiliary 
device since it merely serves to eliminate the tunneling 
contact; the same effect could be achieved by increasing 
the barrier. Our numerical simulations show that /xq and 
fii can be chosen such that the experimentally uncertain 
value of the onset of the pulse to is only of minor impor- 
tance. We determine the optimal values for /io and /ii 
numerically; however, slight deviations from the values 
given below will lead only to minor deteriorations of the 
quality of transfer. 

For the experimentally realistic parameter combina- 
tion Nk — 1.6 f2 we use the Shampine-Gordon rou- 
tine to numerically solve the Schrodinger equation. 
The only difference between the solutions of the N- 
particle Schrodinger equation and its mean field analog 
(the non-linear Schrodinger equation) is a collapse of the 
well-to-well oscillations on the iV-particle level for large 
times which will be followed by a revival |2^; the N- 
particle dynamics displays, on short times scales, a self- 
trapping which is at least as good as on the mean field 
level. We get good results for the transfer if we choose 
/ii = —2.05^ and /io = fJ.i, as shown in Fig. |21 How- 
ever, neither trapping nor transfer is perfect then. For 
a very weakly coupled double well system, we can apply 
much larger constant forces. This will increase the effec- 
tive interaction, and therefore the trapping both before 
and after the tunnehng process (cf. Eq. ((Tn|l ). Choosing 
fj,i = — l.Sfl and fj,o = W ^J^l leads to a transfer of more 
than 99% of the particles (see Fig. [^J. In essence, the 
transfer is achieved by choosing the force-function such 
that the interaction is effectively switched off for the dom- 
inant part of the multi-particle wave- function, according 
to a neutralization of the phase factors © for the most 
important indices v. 

IV. GENERATION OF MESOSCOPIC 
ENTANGLEMENT 

In what follows we exploit this possibility of effectively 
switching off the interaction to produce mesoscopic su- 
perpositions (for various measures of entanglement see 
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FIG. 2: Population transfer as a function of time obtained for 
pulses given by Eq. dJ. Dashed line: numerical solution of 
the equation of motion in the mean field approximation for 
fii — —2.05 51 and fio — ni. Solid line: numerical solution of 
the mean field equation for fii = —1.8Q, and /io = 10 /ii. In 
both cases, the parameter to was chosen such that we obtain 
a worst case-scenario as far as trapping of the final state is 
concerned (Slto = 0.92 and ^Ito = 2.4, respectively). Solutions 
of the A''-particle Schrodinger equation for A'' = 1000 (dotted 
lines) confirm the validity of the mean field approximation. 
Inset: For the second scenario, the maximum and minimum 
value of trapping are plotted as a function of the onset of 
the pulse to (obtained from the mean field equation). The 
transfer to well 2 varies between better than 99% and 99.9%. 



e.g. Ref. |23| and references therein). Our strategy con- 
sists of two steps: First, an initial state has to split into 
two parts. Secondly, we separate both parts by an ap- 
propriate pulse. The initial splitting is already achieved 
by starting with a state which is approximately a Fock 
state \n) with N — n particles in well 1 and n particles 
in well 2. For n — N/2 and no applied force, the prob- 
ability distribution \a^\'^ then evolves automatically into 
a bimodal distribution {cf. e.g. Ref. ^o that the 

state evolving from a Fock state already displays entan- 
glement. A measure for the "quality" of entanglement is 
the distance Av between the two maxima. Fortunately, 
the numerical results of entanglement generation are ro- 
bust against small uncertainties in the preparation of the 
initial state. 

For a superposition of distinct Fock states, the de- 
gree of entanglement does not depend on the phases with 
which the individual states contribute to the total wave- 
function. To characterize such an entangled superposi- 
tion more thoroughly, we perform a gedanken experiment 
by considering the outcome of all possible measurements 
of the number £ of particles found in well 1. We then 
take the average of all measurements with £ < N/2 and 
call the resulting quantity {£) < . Analogously, let {£) > be 
the average over all measurements with £ > N/2. An 
ideal entangled state (|0) -I- \N))/\^ would thus be char- 
acterized by (^)< =0 and (£)> = N. Following the evo- 
lution of an initial Fock state numerically, we find that 
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FIG. 3: Mesoscopic superposition evolving from an initial 
Fock state with 75% of the particles in well 1 and 25% in 
well 2, under the pulse specified by Eq. 11311 . Parameters are 
Nk = 1.6 n and N = 1000. In this plot, probabilities \a^{t)f 
larger than 0.2% are displayed as a function of both time t 
and Fock state index i/. 

{£)y — (£)< oscillates in time. For the particular initial 
state \N/2) and k = 0, we obtain (€)> - (£)< < 0.637 7V. 
Thus, even in the non-interacting case and without any 
appHed force, simply letting an initial Fock state evolve 
will never lead to a perfect superposition of \N) and |0). 
Somewhat surprisingly, for an interacting system in the 
self-trapping regime, better results can be obtained if we 
start with, say, 75% of the particles in well 1 and 25% 
in well 2. Initially, roughly half of the dominant squared 
amplitudes |a^p describing the many-particle wave func- 
tion will move towards lower v, the other half towards 
larger v. To separate these two parts, we apply a pulse 
such that the interaction is effectively switched off for 
that part which shifts towards higher i^; when the other 
half of the wave-function is "reflected" at i/ = 0, the force 
should be such that trapping sets in soon. This intuitive 
scenario is implemented through the force function 

, . _ ( -NKCOs[n{t + ti)] : 0<t<TT/n-ti 

i^Jy^)-\ ■■ t>TT/n-ti 

(13) 

with ^(1 — cos(f2ti)) ~ j^. The constant offset /io is 
again meant to increase the trapping. However, now we 
require strong trapping both for v N and i/ « 0. Thus, 
a parameter /io of the order of Nk will not do — ei- 
ther an offset with a modulus of, say, lONn can be ap- 
plied, or no offset at all. FigureOldepicts the evolution of 
the probabilities |ai/(t)p for finding ly particles in well 2 
under the influence of the force (|I3|I with Nk = 1.6 il 
and iV = 1000, and confirms in detail the above quali- 
tative reasoning: The distribution develops pronounced 
maxima at both low and high v. For the final state, we 
have {£)> - {i)<> 0.73 N . 

Thus, while the final state still is no perfect entan- 
gled state, we have demonstrated that our strategy is 
capable of producing states which at least come close 
to the theoretical ideal. From the numerical estimates 
stated above for ^'^Na, we infer that the production time 



3/n ~ 40 ms is shorter than typical experimental coher- 
ence times 0, 0, . Creating and studying mesoscopic 
superpositions of condensates might greatly improve our 
knowledge of the achievable degree of coherence in meso- 
scopic systems. In a laboratory experiment, the required 
initial Fock state could be prepared by involving a con- 
stant potential offset between the two wells, in order to 
prevent untimely onset of tunneling [Tol |. Our numerical 
data confirm that the resulting state depends only weakly 
on the particular moment when this offset is switched off. 



V. CONCLUSION 

To conclude, we have shown numerically that a suit- 
ably designed time-dependent modulation of a double 
well potential can partially eliminate the effect of the 
interaction between the particles of a trapped Bose- 
Einstein condensate. This finding can be exploited for 
achieving a controlled tunneling transfer of the conden- 
sate from one well to the other, or for preparing meso- 
scopic entangled states. While here we have guessed the 
shape of the modulation by mere inspection of the phase 
factors ||SJ), in order to establish the principal feasibil- 
ity of our strategy, the result might be optimized if an 
analysis of the flnal state is employed, by means of a feed- 
back loop, to design more efficient pulses. Thus, we sug- 
gest to employ the techniques which have already been 
developed and successfully applied to the coherent con- 
trol of molecular dynamics by specifically designed laser 
pulses 20, 21] to the coherent control of condensates in 
specifically modulated double-well potentials. 
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APPENDIX A: VALIDITY OF THE TWO-MODE 
APPROXIMATION 

In order to demonstrate the validity of the two-mode 
approximation, we repeat the mean field calculations for 
a non-linear Schrodinger equation also in a four-mode 
approximation (to do the A^-particle calculations in a 
four-mode approximation for 1000 particles exceeds the 
capacity of most available computers). Throughout the 
appendix, we assume that the individual wells can be de- 
scribed by harmonic oscillators. The validity of the two- 
mode approximation is not restricted to this situation. 
However, to go beyond the two-mode approximation, as- 
sumptions about the potential have to be made. 
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Starting from the non-linear Schrodinger equation 
at \ 2m / 



Ng\^{T,t)\H{r,t) 



with 
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m 



(Al) 
(A2) 



we make the ansatz 

*(r,t) = e"'^«*/''[ci(t)ui(r) + C2(t)ti2(r)] (A3) 



without a time-dependent potential difference {f{t) = 0) 
and without interaction (g = 0), ui^2 are the single par- 
ticle wave-functions for particles localized in the ground 
state of well one/two; wi and W2 are the wave- functions 
of the first excited states. For the weak coupling in 
our physical situation, the solutions of the single-particle 
Schrodinger equation without forcing are given by 

U±{f) ~ [Ul{f) ± U2{f)] 



V2 
1 



Wl(7^ ± W2if)] 



with eigenenergies Eq for u± and for w± , and corre- 
sponding tunneling splittings: 

= E(^ Eq 

Mil = E+ - . 

After inserting Eq. ljA3|l into Eq. ljAl|l and projecting 
the result onto the modes, tedious but straightforward 
calculations for wells which are given as harmonic oscil- 
lators result in the non-linear Schrodinger equation in the 
four-mode approximation: 

-^C2-f Mioe-'^Vm (A4) 
[2ArK|ci|2-)-47VKioMi|' + /i/W] ci 



ICi 



1C2 



irfi — 



[2iV«;|c2p + 47VKio|d2p-Ai/W] C2 

"^d2+/iioe-''"*/(i)ci 
+ [27VAti|di|2-h4iV«;io|ci|2 + ;i/(t)]di 

id2 = -^di-Mioe-'"*/Wc2 

-t- [27VAti|d2p + 4iVAtio|c2p-/i/W]rf2 . 

For the interaction parameters hn = ^ J |mi_2|^ d'^r, 
hni = 1/ |wi_2|''d^r and Hkw = f / |wi,2nwi,2p d^r 
we obtain: 
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FIG. 4: Tunneling beyond the two-mode approximation: frac- 
tion of particles which can be found in the first two excited 
energy levels. The data was obtained by numerically solving 
the non-linear Schrodinger equation in four-mode approxima- 
tion for the same situation as considered before in Fig|5] Top: 
fii — — l.Sr^ and jj,o = 10 ^i. Bottom: /ii = — 2.05f2 and 
Mo = Ml • In both cases the parameters discussed in Sec. ^ 
are used. 



Furthermore, we have: 



(where is the ratio of the oscillator length and half 
the distance of the two potential mimina) . The tunneling 
splittings are given by (c/. 113, 01) ■ 



hill = Ml 



hiUy/2/^d/£,sceM-0-^dyeU , 



^2 



In the derivation of Eq. ljA4|l we used: 



J Ui(f)rMi(r) 


= d 


J U2{r)ru2{r) 


= -d 


j Wi{r)fwi{f) 


= d 


/ W2{r)rw2{r) 


= -d 



where 2c? is the vector connecting the two potential mim- 
ina. 

Within the two-mode approximation, coupling to the 
higher energy levels is neglected and the equations sim- 
plify to: 

ici=~^C2 + 2NK\ci\'^ci+fifit)ci (A5) 

iC2 = -^ci -I- 2iVK|c2pC2 - /U/(i)c2 



6 



T 1 1 r 




FIG. 5: Transfer of the condensate between the weUs: two- 
mode approximation (solid line) versus four-mode approxi- 
mation (dotted line) for fii = — 1.8f2 and fio — 10 ^ii (with 
the same time-dependent potential variation as in Fig.HJ and 
the parameters discussed in Sec. [H] Although the fraction 
of particles transfered to the excited states is of the order of 
one percent (see top curve in Fig. 2J , the quality of the trans- 
fer between the wells differs only very slightly between both 
models (see inset). 



In Fig. 0] the fraction of particles excited to the two 
upper energy levels is plotted as a function of time, for the 
same situation as considered before in Fig. |2 While for 
the case with an initial offset of /ii = —2.05 57 and /io = 
fii [cj. Eq. (|12|l ). the degree of excitation lies below 0.1%, 
it is a order of magnitude larger for /xi = —1.817 and 
/io = 10 /ii. Still, even in this case the transfer between 
the wells differs only very slightly between both models 
(see Fig. 01. 
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